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ABSTRACT 

We present detailed observations of MRC 0116+111, revealing a luminous, mini radio- 
halo of ~ 240 kpc diameter located at the centre of a cluster of galaxies at redshift z = 
0.131. Our optical and multi- wavelength GMRT and VLA radio observations reveal a 
highly unusual radio source: showing a pair of giant (~100 kpc diameter) bubble-like 
diffuse structures, that are about three times larger than the analogous extended radio 
emission observed in M87 - the dominant central radio galaxy in the Virgo Cluster. 
However, in MRC 0116+111 we do not detect any ongoing Active Galactic Nucleus 
(AGN) activity, such as a compact core or active radio jets feeding the plasma bub- 
bles. The radio emitting relativistic particles and magnetic fields were probably seeded 
in the past by a pair of radio-jets originating in the AGN of the central cD galaxy. 
The extremely steep high-frequency radio spectrum of the north-western bubble, lo- 
cated ^100 kpc from cluster centre, indicates radiation losses, possibly because having 
detached, it is rising buoyantly and moving away into the putative hot intra-cluster 
medium. The other bubble, closer to the cluster centre, shows signs of ongoing particle 
re-acceleration. We estimate that the radio jets which inflated these two bubbles might 
have also fed enough energy into the intra-cluster medium to create an enormous sys- 
tem of cavities and shock fronts, and to drive a massive outflow from the AGN, which 
could counter-balance and even quench a cooling flow. Therefore, this source presents 
an excellent opportunity to understand the energetics and the dynamical evolution of 
radio-jet inflated plasma bubbles in the hot cluster atmosphere. 

Key words: galaxies: clusters: individual: MRC 0116+111 - galaxies: active - galax- 
ies: acceleration of particles magnetic fields - X-rays: galaxies: clusters - radio contin- 
uum: general. 



1 INTRODUCTION 

Multi-wavelength observations of galaxy clusters have pro- 
vided substantial evidences for supermassive black holes 
(SMBH) residing in the nuclei of giant elliptical galaxies 
which are found near the cluster centres and dense galac- 
tic environments. The extremely large mass of these black 
holes (Mth ~ 10 6 — W 9 Mq) raises several questions: what 
processes regulate the growth of these SMBH? why and 
how only some of them turn-on their radio jets? what is 
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the source of fuel for the central-engine? and what environ- 
mental impact could a massive black hole and its energetic 
outflows may have on the host galaxy and the surround- 
ing gaseous intra-cluster medium (ICM)? Some important 
clues to these questions are provided by radio and X-ray ob- 
servations of galaxy clusters which show that majority (as 
many as ~70 - 100 percent) of these central giant ellipti- 
cal galaxies (cD galaxies) may become radio-loud - ejecting 
plasma and magnetic fields in the form of powerful radio- 
j ets - when t h ey are located at the focus of a 'cooling-flow ' 
i|Burnsl Il99d : iBagchi fc Kapahil Il994 iMittal et all I2OO8I ). 
The probability of finding a central radio source increases 
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in dir ect proportion to decrease in cooling time of the inner 
ICM l|Bagchi fe Kapahilll994l ; lMittal et alll2008l ). A cooling- 
flow is established when the central radiative cooling time 
(tcooi) of the ICM of a galaxy cluster is quite short (t coo i 
<1 Gyr << Hq , the Hubble time). In the absence of 
a central heating mechanism, the pressure balance is dis- 
turbed and to restore the hydrostatic equilibrium, a steady 
convergent flow of cooling gas towards the cluster centre is 
setup, in which the gas cools below the X-ray temperature 
and accretes onto the cent ral elliptical galaxy (|Fabianll 19941 ; 
iPeterson and Fabianl [20061 '). 

Strongest spectral signatures of this cooling gas are 
low energy X-ray lines below 1 keV from various ioniza- 
tion states of Iron (mainly Fe L complex). Early Einstein 
observations of nearby radio-loud cD galaxies in cooling- 
flow clusters such as M87 (Virgo-A) and NGC 1275 (3C84, 
Perseus-A) did find th e expected Fe L lines, thus support- 
ing the cooling model l|Canizares et al.| 1982). Evidence for 
cooling gas was also seen in other wavelengths, although at 
rates lower than expected; H a e mitting filaments of ion- 
ized, warm gas dHeckmanl [l9 81), star formation activity 
(McNamara and Q'Conellll989l) and radio observa tions of 
CO molecular lines (|Edgell200ll ; ISalome et allbOOSl ). 

In recent years more detailed X-ray observations with 
Chandra and XMM-Newton telescopes have revealed a sur- 
prising and somewhat puzzling aspect of cooling flows; they 
showed far less cooling below X-ray emitting temperatures 
than expected, altering the previously accepted p i cture of 
cooling flows (cf. review by IPeterson and Fabianl i|2006h ). 
Unless gas is thermally supported, radiative cooling leads 
to a 'cooling catastrophe', i.e., an inexorable inflow of cold 
gas onto the central galaxy. To prevent this, some heating 
mechanism was required to raise gas temperature above ~ 2 
keV, suppressing the cooling flow. Although several such 
mechanisms were discussed, the most promising heating pro- 
cess turned out to be the energy injected into the ICM by 
powerful radio jets emanating from AGNs in central galax- 
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ies of clusters and groups (e . g „. iBinnev and Taborl 
IChurazov et all (1200 il. 120021); iBruggen and Kaiserl 
and see iMcN amara and Nulsenl (|2007l) for a review). This 
AGN-heating mechanism is also found to be particularly ef- 
fective in suppressing the star formation activity, that one 
expects to be prolific around the brightest cluster galaxies 
(BCGs) due to their preferential location in clusters. These 
BCGs should be accompanied by intense star formation ac- 
tivity and have blue colours, neither of which is observed 
and energy input from AGN outbursts in clusters may be 
especially needed to explain the observations. 

One of the canonical models of AGN posits accre- 
tion of gaseous medium as fuel for the nuclear black 
hole, such that AGN outflows are powered by gravita- 
tional binding energy released by the infalling matter 
|Begelman. Blandford and Reesll 19841 ). The observed strong 
association of an AGN in the central galaxy and the sur- 
rounding cooling-core lends good support to this model, 
which suggests that the black hole 'cent ral engine' is pos- 
sibly fueled by accretion of cooling gas (IBagchi fe Kapahl 
1 19941 ; lAllen et all 120061 ; iMittal et all 120081). with the flow 
rate self-regulated by AGN -heating |Churazov et alj l200ll ; 
ICattaneo and Tevssier|[2007l ) , indicating a complex feedback 
loop with tight coupling between a central black hole and 
the surrounding gas of cooling core. Although this model is 



plausible and widely used, many details of how this feedback 
process works remain far from clear. 

When the radio jets emerging from the central black 
hole (AGN) interact with the dense thermal plasma of the 
ICM, two bubble-like lobes of non-thermal plasma are in- 
flated, which are filled with relativistic particles and mag- 
netic field and thus become visible in radio observations. 
Such bubbles were first proposed by iGull and Northoverl 
(|l973l ) and later identified in sev eral clusters using ra- 
dio a nd X-ray observations (see iMcNamara and Nulsenl 

(2007) ). The most clear example of this phenomenon are the 
non-thermal bubbles in clusters MS0735.6+7421, Hydra-A, 
Abell 2052, Perseus and others, showing an unusually large 
and energetic pair of r adio e mitting, X-ray dark c avitie s 
(e.g.. IMcNamara et ail (|2000| , IgOoBT) ; Bla nton et alj (|2003h ; 
IMcNamara and Nulsenl (|2007l ); 1wTse et all (120071 ))? Episodic 
(on-off) activity of radio jets injects non-thermal radio bub- 
bles which may heat the ICM via weak shocks, and addition- 
ally these plasma bubbles are responsible for the mechanical 
(PdV) work done on the ICM for heating it, which is one of 
the favored mechanism of AGN-ICM feedback. 

Several examples of these bubble and cavities, with di- 
ameter ranging from 1-100 kpc, have been found in the 
hot a tmosphere su r round i ng galaxies in groups and clusters 
(e.g- lBirzan et all (|2004l ); lDunn and Fabianl j|2004l )). In few 
sources the energy involved is ~ 1 60 ~ 62 erg, the mo s t pow - 
erful radio outbursts known ( e.g., IMcN amara et alj (2005); 
IMcNamara and Nulsenl (|2007l )). The energy involved is large 
enough to strongly affect or even quench any cooling flow, 
and to drive large-scale outflows that redistribute and heat 
the gas on cluster-wide scales. The radio lobe plasma fill the 
cavities, which shows that expanding jets have displaced the 
ICM thermal gas, excavating X-ray dark cavities and pos- 
sibly shock h e ating t he surrounding ga s eous mediu m (e.g., 
iFabian et all (120031 ); IMcNamara et all (|2000l . 120051 )). This 
heat input into the ICM may be large enough to counter- 
balance the cooling loss, and may even stop the accretion 
flow of matter onto the SMBH, thus starving the central en- 
gine of its fuel. Therefore, bubble-like diffuse radio sources 
residing near cluster centres are extremely effective probes 
of the poorly understood physics of the radio galaxy-ICM 
feedback process. 

Also directly related to these issues is the detection of 
diffuse non-thermal synchrotron radio plasma and magnetic 
fields mixed with the ICM of clusters. Approximately 30 
galaxy clusters show giant radio haloes - centrally located 
synchrotron radiation extended on mega-parsec scales with 
a regular morphology, similar to that observed in X-rays. 
A few galaxy clusters also show peripherally located radio 
relics which have elongated morphology and that may ex- 
tend over 100 kpc - 1 Mpc scales. The radio emission from 
giant radio haloes and radio relics is not directly connected 
to the galaxies in clusters and both are mainly observed in 
merging clusters, which indicates that radio emitting rela- 
tivistic electrons may have been accelerated in- situ in merger 
induce d shocks and/or turbulence (cf. review iFerrari et all 

(2008) and references therein). In the case of relics, an al- 
ternative model has been proposed, which invokes adiabatic 
compression of relic radio lobes by cluster merger shocks 
j|Ensslin and Gopal-Krishnall200ll ). 

A third remarkable category of large scale diffuse radio 
sources are the radio mini-haloes. Unlike the mega-parsec 
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scale giant radio haloes and relics (which are found to be as- 
sociated with cluster merger phenomenon and hence rarely 
found in cooling core clusters), radio mini-haloes are not 
only a few times smaller in size (~100 - 500 kpc, comparable 
to size of the cooling core) but they are seen to surround the 
centra l radio galaxies of some cooling core clusters (see r e- 
views: IFeretti and Giovanninil (120081 ); iFerrari et all i|2008l )). 
Radio mini-haloes are rare objects, have low surface bright- 
ness and a steep spectral index (a^-10 The fact that 
mini-haloes are mostly observed at the centres of cooling- 
core clusters indicates that their origin and evolution are 
connected to the energy feedback from the central AGN 
into the ICM via radio jets and cooling/heating processes. 
Thus, while mini-haloes share with gaint-haloes and relics 
the properties of a steep radio spectrum and very low surface 
brightness, their physical origins are probably different. This 
is corroborated by the finding that unlike the mini-haloes, gi- 
ant haloes have remarkably similar synchrotro n emissivities, 
albei t the two are similar morphologically (|Murgia et al.l 
120091 ). Radio mini-haloes are still a very poorly understood 
phenomenon and in observational terms, they being roughly 
centred at the dominant cluster radio galaxy makes their de- 
tection particularly challenging. Three early examples of ra- 
dio mini-ha loes are assoc iated with t he cooling flow clusters: 
Perseus- A l|Pedlarl H99fjh . A 2142 (iGiovannini and Ferettil 
|2000| ) and RX J1347.5-1145 (|Gitti et al.ll2007l ). Continued 
sear ches have so f ar revealed no more than 10 mini-h aloes 
fsee lFerrari et all (120081 ) ; IFeretti and Giovanninil (|2008h and 
references therein) . Due to their very steep spectrum, radio- 
haloes in clusters are difficult to detect at higher frequen- 
cies, and therefore a low frequency radio tel escope such as 
the GMRT is well suited for their study (e.g.. lBrunetti et al.l 
i|2008MVenturi et alj i|2008l )). 

The focus of present paper is on a diffuse halo-like radio 
source MRC 0116+111 and its radio and optical properties, 
as revealed in our detailed observations, their astrophysi- 
cal implications and impact on the above topics. Earlier a 
brief report on this radio source matching the characteris- 
tics of a mini-halo, was present ed by some of us, based o n 
VLA and GMRT observations (jGopal-Krishna et al-lliool ). 
These radio observations revealed an amorphous source of 
an ultra-steep radio spectrum (a ~ —1.3) and size ~ 1.5'. 
This Molonglo Reference Catalogue source MRC 0116+111, 
had been resolved on arc-minu te scale in the Ooty L unar 
Occultation Survey at 327 MHz ijjoshi and Si ngal 1980]). An 
R-band image taken with the ESO New Technology Tele- 
scope (NTT) showed that this amorphous radio source is 
centred at the dominant member of a galaxy group having 
at least 3-4 members for whic h NTT/EMMI spectroscop y 
yielded redshifts close to 0.131 (|Gopal-Krishna et al.ll2002l ). 
More recently, this system has been listed in a catalogue of 
distant clusters of galaxies (jLopes et alJ 1200*41 ). In order to 
clarify the nature of this peculiar radio source and its re- 
lation to the brightest cluster galaxies, we have imaged it 
freshly with GMRT at 240 MHz, 621 MHz and 1.28 GHz 
with high sensitivity and resolution, and have also carried 
out an improved analysis of the existing VLA snapshot data 
at 1.4 and 4.8 GHz. Additionally, we have taken its deep 
B, V,R,I CCD images with the 2-mt optical telescope lo- 
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Figure 1. Upper panel: GMRT 1.28 GHz map of 
MRC 0116+111 shown with contours (levels: ±0.24, 0.48, 0.96, 
2, 4, 8 mjy/beam, noise r.m.s. 0.08 mjy/beam, beam: 5" FWHM 
circular, plotted inside box) overlayed on IGO R-band image. No 
AGN (radio core) is visible down to ~ 1 mjy/beam flux den- 
sity limit, and no radio jets or lobes are detectable either. Lower 
panel: GMRT 240 MHz map shown with contours (levels: ±4, 
8, 16, 32, 64, 128 mjy/beam, noise r.m.s. 1.35 mjy/beam, beam: 
12.15" X 10.31" FWHM at 65.4° PA, plotted inside box). In the 
background is shown the IGO R-band image of the galaxy cluster. 



cated at the IUCAA Girawali Observatory (IGO), Pune, 
India. 

In this work we have used a Hubble constant Ho = 
70 km/(s Mpc), and A-CDM 'concordance' cosmology with 
Am = 0.27 and J1a = 0.73, which results in a luminosity 
distance Dl = 617.2 Mpc, and linear scale of 2.34 kpc per 
arcsec for a redshift z = 0.131 of MRC 0116+111. 
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2 OBSERVATIONS AND DATA ANALYSIS 
2.1 Radio observations 

Multi-frequency Giant Meterwave Radio Telescope 
(GMRT0 observations were carried out on December 
1 and 2, 2007 at 1.28 GHz (A = 23.4 cm), 621 MHz 
(A = 48.3 cm) and 240 MHz (A = 125 cm) using the 
128 channel FX correlator. GMRT is an earth-rotation 
aperture-synthesis ar ray of 30 tele s copes operating in 150 
- 1420 MHz range ijSwarup et all Il99ll ). The data were 
recorded in the spectral-line mode with 128 channels per 
sideband, resulting in a spectral resolution of 125 kHz per 
channel. The 621 and 240 MHz observations were taken 
simultaneously with the co-axial dual-band feed system of 
GMRT. At 240 MHz we used a bandwidth of only 6 MHz 
out of total 16 MHz available in one sideband (USB), and 
at 621 MHz full bandwidth of 16 MHz in one sideband 
(USB) was used. At 621 MHz only the right-handed circular 
polarization (RR) and at 240 MHz only the left-handed 
circular polarization (LL) was recorded. Both right and left 
circular polarization (RR, LL) data were recorded at 1.28 
GHz, and bandwidth of 16 MHz was used in one sideband 
(USB). During each observing run a sequence of ~20 to 
30-minute scans of the source were taken, interspersed 
with the VLA phase calibrator 0204+152 (4C +15.05). 
The telescope gain, bandpass response and flux density 
scale were calibrated by observing the primary flux density 
calibrator 3C 48, and boot-strap ping the flux dens ity scale 
to the standard 'Baars-scale' (|Baars et al.l Il977f ). After 
excising faulty data and RFI, mainly affecting low frequency 
observations, the calibrated data were transformed into the 
image plane using the standard routines in AIPS, including 
a few cycles of self-calibration. By applying suitable tapers 
to the visibilities, final maps with synthesized Gaussian 
beams of FWHM ranging 5 to 6" at 1280 and 621 MHz, 
and ~10" at 240 MHz were obtained (see Fig. [T] and 

Fig. ED . 

The first of Very Large Arra^ observations took place 
on June 22, 1992 in C-band (DnC array configuration, 
4.88 GHz and 4.83 GHz centre frequencies, using two 
IFs of 50 MHz BW each). The total integration time on 
MRC 0116+111 was about 10 minutes and VLA calibrator 
B0235+164 was used for phase calibration. VLA calibra- 
tor 3C 48 was also observed for amplitude calibration. The 
second VLA observation was made on June 12, 1997 in L- 
band (CnB configuration, 1.46 GHz and 1.38 GHz centre 
frequencies, using two IFs of 50 MHz BW each). The total 
on-source integration time was about 38 minutes and VLA 
flux calibrator 3C 48 was used for both amplitude and phase 
calibrations. For data reduction we used standard routines 
available in AIPS software. Further editing, and imaging 
with a few cycles of self-calibration resulted in good quality 
radio maps. After deconvolution, the final maps were re- 
stored with synthesized Gaussian beams of FWHM 12" at 
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National Science Foundation, operated under cooperative agree- 
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both the frequencies (see Fig. [2}. Table 1 gives the details 
of observing runs with the different radio telescopes and the 
sensitivities and resolutions that were achieved at various 
observing frequencies. 

For obtaining a spatially resolved spectral index in- 
formation, we combined GMRT 240 and 621 MHz images 
for a low frequency spectral index map, and VLA 1.4 and 
4.8 GHz images for a high frequency spectral index map. 
These are shown in Fig. HJ] For this purpose, radio data of 
MRC 0116+111 were specially processed (if required, by giv- 
ing suitable tapers in the visibility plane) so as to obtain a 
uniform angular resolution of 12" (FWHM beam size) at all 
the frequencies. When combining the radio images we have 
included only the pixels with amplitudes ~ 3.5 times above 
the r.m.s. noise, to prevent spurious structures from appear- 
ing in the spectral index maps. These cut-offs were given at 
values 4 mjy/beam, 0.5 mjy/beam, 0.27 mjy/beam and 
0.15 mjy/beam on the 240, 621, 1425 and 4860 MHz maps 
respectively. 

2.2 Optical imaging and spectroscopic 
observations 

The diffuse radio-halo source MRC 0116+111 was ear- 
lier identified by us to be centred on the dominant 
(cD-like) galaxy of a po or group with 3-4 members 
l|Gopal-Krishna et al.l I2002T ). Recently, at the position 
of this radio source a relatively p oor galaxy cluster 
NSCS J011904+112133 was identified bv lLopes et ail (|2004h 
in their survey of clusters from DPOSS-II photographic 
plates, with a photometric redshift estimate 0.16. However, 
From our NTT sp ectroscopy we placed t he clu ster at a red- 
shift of 0.131 (see iGopal-Krishna et all (|2002h and below). 
Other than this no other optical or spectroscopic data was 
found in the literature. We carried out optical broadband 
(B, V,R,I) imaging observations of the host cluster using the 
recently commissioned 2-mt telescope at the IUCAA Gi- 
rawali Observatory (IGO), located about 80 km from Pune, 
India. The observations were taken on the night of 31st De- 
cember 2007 using the IUCAA Faint Object Spectrograph 
and Camera (IFOSC) mounted at the Cassegrain focus. 
IFOSC employs an EEV 2Kx2K, thinned, back-illuminated 
CCD with 13.5 pirn pixels. The spatial sampling scale at the 
detector is 44/xm arcsec -1 giving a field of view of about 
10.5 arcminutes on the side. During observation the seeing 
conditions were good, but with fair to occasionally poor sky 
transparency. 

We took several frames in each filter with total exposure 
time sufficient to achieve adequate photometric accuracy. 
Photometric zero points and colour transformations were 
defined by ob serving Landolt's standard field Rubin 149 
i|Landoltlll992l ). The analysis of the CCD frames was done in 
a standard way using the IRAF software. The /-band frames 
suffered from strong fringing, which could not be completely 
removed by flat fielding; therefore we report only the B, V 
and R band data in this paper. A composite B,V,R colour 
image of the host galaxy cluster is shown in Fig. [3] Several 
elliptical galaxies are visible within the contours of GMRT 
and VLA radio maps (Fig. [T] and [2} with a giant cD-like 
galaxy of m v = 17.69 mag (M v — —21.38 mag) seen in 
the centre at the position RA = 01 h 19 m 1.69 s , Dec - +11° 
23 18.4" (J2000). 
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Table 1. Details of GMRT and VLA observations of MRC 0116+111 



Frequency (MHz) 



210 



621 



1280 



1425 



4860 



Telescope 

Date of Observation 
Flux Calibrator 
Bandwidth (MHz) 
Integration time per visibility 
FWHM synthesized beam 
RMS map noise (mjy/beam) 
Integrated flux density (Jy) 



GMRT 
1 Dec. 2007 
3C 48 
6 

16.9 s 
12.15" x 10.31" 
1.35 
1.15 ± 0.05 



GMRT 
1 Dec. 2007 
3C 48 
16 
16.9 s 
6" circular 
0.15 
0.47 ±0.01 



GMRT 
2 Dec. 2007 
3C 48 
16 
16.9 s 
5" circular 
0.08 
0.19 ±0.01 



VLA 
12 June 1997 
3C 48 
50 
10 s 
12" circular 
0.08 

0.140 ±0.005 



VLA 
22 June 1992 
3C 48 
50 
10 s 

12" circular 
0.04 

0.035 ± 0.005 



Table 2. Observed and calculated parameters 



Spectroscopic redshift (central cD galaxy) 
Apparent magnitudes t (central cD galaxy) 
621 MHz Radio Luminosity (1/621 MHz) 
1400 MHz Radio Luminosity (Li.4GHz) 
Bolometric Radio Luminosity (L ra dio) 
(over 10 MHz - 10 GHz range) 
'Break' frequency in spectrum 
Elliptical North-west bubble dimension 
(projected semi-major and minor axes) 
Circular South-east bubble dimension 
(projected circular radius) 



0.1316 

m B = 18.79, m v = 17.69, m R = 16.69 
1.21 x 10 25 W Hz" 1 
4.57 x 10 24 W Hz" 1 
3.64 x 10 34 W 

~ 400 MHz 

~27.5" X17.5" (= 64.4 kpc x41 kpc) 
~35" (= 82 kpc) 



f Corrected for galactic extinction (NASA Extragalactic Database). 



For spectroscopy we used the ESO/La Silla 3.6-meter 
New Technology Telescope (NTT) and EMMI (ESO Multi- 
Mode Instrument). A 2x10 minutes low resolution slit- 
spectrum (slit width 1.5" , spectral resolution 3.7A) was 
taken with the grism-3 optics on the EMMI (October 1991) 
which gave a redshift z = 0.1316 for the brightest central cD 
galaxy (Fig. [4j) , based on the absorption lines of Na (A5893), 
Mgb (A5169), H/3 (A4861), G-band (A4304), and the Ca H,K 
break. A probable weak [Oil] (A3727) emission line is visible 
at the left-most edge of the spectrum, where the S /N ratio is 
comparatively lower. The second brightest elliptical galaxy 
~ 15" south of cD (Fig.[3| was found to have a close redshift 
2 = 0.1309, suggesting a group-like environment around the 
central cD. 



3 PHYSICAL PICTURE OF THE RADIO 
SOURCE 

3.1 A twin bubble-like diffuse radio-halo 

morphology: A 'quiescent' radio galaxy? 

MRC 0116+111 shows a highly diffuse 'halo' or 'bubble' like 
structure which bears close resemblance to the amorphous 
structure of so-call ed 'radio mini-haloes' (see Sect. 1 and 
iFerrari et all ((2008)). The physical size, radio spectral index 
and morphology of MRC 0116+111 are broadly consistent 
with the properties of a cluster radio mini-halo. In Table. 3 
we compare the physical properties of MRC 0116+111 with 
a few other well known radio mini-haloes, including the well 
known radio galaxy Virgo-A (M87). We note that the ra- 
dio luminosity of MRC 0116+111 at 1.4 GHz, Li. 4 ghz = 
4.57 x 10 24 W Hz -1 , is quite large, comparable to the lumi- 



nous radio-haloes in Perseus and A2390 clusters. Also, its 
bolometric radio luminosity (over 10 MHz - 10 GHz range), 
iradio = 3.64 x 10 34 W, close to the divide between FRI and 
FRII class radio sources, would place it amongst the most 
luminous radio haloes known. 

GMRT 621 MHz radio map of MRC 0116+111 (Fig. EJ 
and its other radio maps (Fig[T] [5J clearly show that on 
larger scales (~100 kpc) the radio-halo emission comprises 
of two diffuse 'lobe'-like structures surrounding the central 
group of galaxies, which we suggest are actually a pair of 
giant, radio emitting plasma bubbles, filled with relativistic 
particles (cosmic rays) and magnetic field. In Figure [7] the 
dotted lines delineate this pair of radio bubbles, having the 
projected dimensions ~27.5" X17.5" (= 64.4 kpc x41 kpc) 
- the semi-major and minor axes of ellipse for the smaller 
elliptical north-west bubble, and ~35" (= 82 kpc) radius for 
the larger spherical bubble to south-east of the cluster cen- 
tre. 

The total spatial extent of radio emission is ~240 kpc 
which apparently surrounds the dominant brightest cD-like 
elliptical galaxy of the cluster (Fig. |3J . However, based on 
the detection sensitivities achieved in our present GMRT 
and VLA observations (see Table 1), there is no strong evi- 
dence that the extended radio emission emanates from this 
or any other galaxy enclosed within the radio contours. The 
highest resolution GMRT 1.28 GHz map (Fig. [TJ» and other 
radio maps have not detected a compact active radio nu- 
cleus (AGN) and radio-jets directly feeding into the 100-kpc 
scale twin plasma lobes/bubbles of MRC 0116+111. In this 
respect MRC 0116+111 is very unusual and it differs from 
other well known mini haloes, including canonical M87 and 
Perseus- A (jFerrari et al.ll2008l ; [c assano. Gitti and Brunetta 
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RISHT ASCENSION <!J239e> 

Figure 2. VLA radio maps of MRC 0116+111 overlayed on 
IGO R-band optical image of the galaxy cluster. The up- 
per panel shows L-band (1425 MHz) map drawn with con- 
tour levels ±0.24,0.48,1,2,4,8 and 12 mjy/beam and lower panel 
shows the C-band (4860 MHz) map drawn with contour levels 
±0.2,0.4,0.8,1.6,2.2 and 3 mjy/beam. The resolution is 12" for 
both the images, shown by the FWHM of the beam plotted in 
the lower left corner. 



|200S| ; lOwen. Eilek and Kassiml l200d ; iPedlad Il990l ). which 
have detectable powerful AGN cores and jets. The high reso- 
lution GMRT 1.28 GHz image of MRC 0116+111 also do not 
show any other compact AGN-related background sources 
blended with the radio-halo emission, which clearly rules 
out its origin from such a contamination (Fig. [TJ). 

Therefore, this suggests to us a physical picture in 
which the relativistic particles and magnetic fields of 
MRC 0116+111 halo were initially seeded by a pair of radio- 
jets originating in AGN activity in the past, most likely from 
the nuclear SMBH of the central cD galaxy, as revealed by 
the pair of plasma bubbles straddling this galaxy. However, 
at present MRC 0116+111 resembles a quiescent or 'fossil' 




Figure 3. A B, V,R false-colour composite of the central region 
of galaxy cluster hosting MRC 0116+111. The region shown is 
about 8.3' x7' in size and north is on top and east is to the left. 
This image is made using the IUCAA Girawali Observatory (IGO) 
2-meter telescope. 
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Figure 4. EMMI (ESO Multi-Mode Instrument) low resolution 
(spectral resolution 3.7A) slit-spectrum of the central cD galaxy 
in MRC 0116+111 cluster taken with grism-3 optics on the New 
Technology Telescope (NTT). The abscissa shows observed wave- 
lengths and ordinate, the spectral flux on an arbitrary scale. The 
prominent spectral lines are marked along with their rest wave- 
lengths. The spectrum gives a redshift z = 0.1316 for this galaxy. 



source - in the sense that its twin extended lobes/bubbles 
of radio plasma not energized by active jets from an AGN 
- as observed in classical radio galaxies. This suggests that 
the SMBH embedded in the nucleus of the cD galaxy, which 
was once active, probably has turned-off its radio jets some- 
times back in the past. Only a handful of such 'fossil' or 
'dying' radio galaxies are known so far, which are more eas- 
ily detected at low radio frequen c ies due to their very s teep 
spectra (e.g.. |Parma et al] l|2007h ; [jamrozv et al.l(|2004l )). It 
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Integrated Radio Spectrum of MRC 0116+111 
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Figure 5. The integrated radio spectrum of MRC 0116+111 be- 
tween 74 MHz and 4.9 GHz. Piecewise linear least-square fit to 
the data points shows a sudden spectral steepening beyond the 
'break' frequency v\y r ~400 MHz. The low-frequency spectral in- 
dex is o; = — 0.55(±0.05) and the high-frequency spectral index is 
a h = -1.35(±0.06). The GMRT and VLA observed flux density 
values are listed in Table 1 . Other data points are tak en from lit- 

327 MHz 



erature as follows: 74 MHz [V LSSi lCohen et ail J2007 



iGopal-Krishna et alJ ll2002lll. 408 MHz [Molonglo: lLarge et all 
dl98ll )1, 1400 MHz [NVSS:ICondon etHI l ll99ct) l. 2700 MHz [Ef- 
felsberg: [Gopal-Krishna et al l ll2002l) l and 4.8 GHz [GreenBank: 
iBecker. White and Edwards! lll99ll ll. 



is possible that we have caught the central massive black 
hole, which created MRC 0116+111, during a low or even 
quiescent state of its activity cycle. If the heat input into 
the ICM by these expanding radio bubbles is large enough 
to significantly offset the cooling loss, the flow of accreting 
matter onto the SMBH may stop, thus starving the 'central 
engine' of its fuel. So we expect that this radio bubble-fed 
energy exchange should give rise to an episodic triggering of 
the AGN itself, establishing a self-regulated activity pattern 
between cooling flow and AGN activity. We point out that 
this cyclic process, even if in operation, does not imply that 
some form of in situ re-acceleration of particles in shocks or 
turbulence does not take place in the central cooling core 
and the extended lobes of MRC 0116+111 - a scenario, in 
fact supported by both the integrated radio spectrum and 
the spectral index maps discussed below. 



3.2 Possible location of MRC 0116+111 in a 
low- mass poor-cluster /group environment 

According to our present understanding of the origin and 
evolution of mini-haloes, their diffuse, steep-spectrum halo 
or bubble like morpholgy results from strong interaction 
of radio jets with the surrounding dense and cool gaseous 
medium of the cluster core. At present there are no avail- 
able good X-ray data for this system and therefore we lack 
information about the physical state of ICM in the core of 
cluster hosting MRC 0116+111. Such data are required for 
quantifying accurately the energy input rate by the AGN 
jets into the plasma bubbles of mini-halo as well as their 
dynamics and interaction with the ICM. The X-ray lumi- 
nosity is also a good indicator of the optical richness and 
mass of the host system. We find that the mini radio-halo 



MRC 0116+111 is located in an under-exposed region of 
ROSAT All Sky Survey (RASS), which only allows us to de- 
termine an upper limit of <1.4x 10 43 erg s _1 for its X-ray 
luminosity in 0.5 - 2 keV band. 

However, based on the mass-X-ray luminosity re- 
lation, the upper limit on the X-ray luminosity of 
MRC 0116+111 translates to an upper limit for the grav- 
itatio nal mass of its host clu s ter/gr o up of about 10 14 M p, 
(e.g., iReiprich and Bohringeri (|2002l ); IStanek et all l|2006h ). 
This suggests that MRC 0116+111 is possibly located at 
the centre of a poor-cluster or a rich-group environment, 
which is also corroborated by the (ph otographical l y esti - 
mated) low galaxy density reported by lLopes et al.l (|2004T ) 
for the galaxy cluster NSCS J011904+112133, which hosts 
MRC 0116+111. This is a rather surprising result which 
makes MRC 0116+111 special because in general, all other 
known mini radio-haloes are found within X-ray bright cool- 
ing cores (Table 3) with X-ray luminosity >10 44 erg s , 
and some of the most powerful systems may have luminos- 
ity 



-10 45 erg s" 1 . 



3.3 Spectral index plot: The 'break' frequency 
and radiative ageing analysis 

As these 100 kpc scale giant radio-emitting magnetized bub- 
bles in MRC 0116+111 suggest their origin in intermittent 
AGN activity, we can estimate the timescale of such out- 
bursts from the shape of total radio spectrum and the spec- 
tral index maps. Spectral index represents a powerful tool 
to study the properties of the relativistic electrons and the 
magnetic field in which they radiate, and to investigate the 
connection between the electron energy and the ICM. From 
our multi-frequency radio imaging and published data, we 
have derived the integrated spectrum (Fig. [SJ) and spatially 
resolved spectral index maps (Fig. [6}. Low frequency ( ~100 
- 400 MHz) spectra are important to determine the slope of 
electron energy distribution and preserve a record of past 
AGN activity, while the high frequency spectra (~1 - 10 
GHz) give information on the diffusion and ageing of rel- 
ativistic particles, and any particle (re) acceleration mecha- 
nism that may be in operation. 

The integrated radio spectrum between 74 MHz and 
4.9 GHz frequencies shows a strong downward curvature 
(Fig. [5}. Piecewise linear least-squares fit shows an onset 
of spectral steepening beyond the critical (also known as 
'break') frequency v c ~400 MHz. The low-frequency spec- 
tral index is ai — — 0.55(±0.05), while the high-frequency 
spectral index attains a slope aj, = — 1.35(±0.06). The steep 
spectral index of MRC 0116+111 is consistent with the mini 
radio-halo mor phology whic h usua lly have steep radio spec- 
trum (a < -1. iFerrari et alJ {2008)). If the relativistic elec- 
trons are all injected in a single energetic outburst with a 
power-law energy distribution, N(E)dE = NoE~ p dE (here 
iVo is amplitude and N(E)dE the differential number den- 
sity within energy E to E + dE, and p the power law in- 
dex), subject to strong radiation losses in synchrotron and 
inverse Compton processes, the radio spectrum would cut- 
off rapidly to zero for frequencies above a critical frequency 
v > u c , which shifts to lower values with time. In contrast, 
if episodic injection of fresh particles, or re-acceleration of 
existing particles takes place in the source, only a spectral 
steepening by ~0.5 beyond the critical frequency occurs. In 
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MRC 0116+111 beyond the critical frequency v c the radio 
spectrum steepens due to radiative losses, but no rapid cut- 
off is seen (Fig. [5} , which suggests an episodic or multiple 
particle injection process rather than a single energetic out- 
burst event. 

The spectral age t av derived from the synchrotron radio 
spectrum is given by (|Slee et alj|200ll ): 



1.59 x 10 J 



yr, 



(1) 



where the magnetic field B is in fj,G, the critical frequency 
u c in GHz and Bqmb = 3.2 x (1 + z) 2 /iG is the equivalent 
'magnetic field' of the cosmic microwave background (CMB) 
radiation at redshift z. This formula assumes isotropic pitch- 
angle distribution, no expansion losses and a uniform mag- 
netic field which remained unchanged over the radiative age. 
For v c = 400 MHz and magnetic field range B= 1 — 10 
/j,G, the electron spectral age is t sp — (1.33 — 0.64) x 10 s 
y, which can be taken as the elapsed time since last injec- 
tion of relativistic particles in the source. The largest uncer- 
tainty in estimation of spectral age is imposed by magnetic 
field B. Presently we do not have X-ray observation of this 
source, but from the observed mini radio-halo morphology, it 
is likely that MRC 0116+111, in common with other known 
mini radio-haloes, would be located at the centre of a cluster 
cooling-core, which in general have stronger ICM magnetic 
fields of order B= 10 — 30 /xG compared to non cooling-cor e 
clusters (|Govoni and Ferettill2004 ISoker and Sara zin 1990). 
Therefore t sp ~fewxl0 7 y is a more likely (globally aver- 
aged) radiative time-scale for MRC 0116+111. 

The injection of relativistic particles in the medium can 
take place via an episodic on-off mechanism of the radio- 
jets driven by the central AGN (even though no radio core 
or jets are visible down to ~1 mjy/beam level), or via a 
Fermi- type shock/turbulent re-acceleration mechanism. The 
requirement is that cooling time t coo i « t 3p should be longer 
than either the time interval tint between two successive nu- 
clear outburst, or the acceleration time-scale t acc for elec- 



trons emitting at frequencies v 



> , 



i.e., t sp > ti n t or 



t sp > tacc- It should be added that a lack of spectral cut-off 
upto the highest observed frequency v — 4.86 GHz, where 
electron cooling time-scale t coo i « 1.8 x 10 7 y (B = 10 fJ-G), 
puts a severe constraint on the acceleration time-scale, which 
should be shorter than this, irrespective of the acceleration 
mechanism. 



3.4 Spectral index maps: particle re-acceleration 
and a buoyantly rising plasma bubble? 

Here we present the spatially resolved spectral index maps 
which reveal interesting details of the radiative energy loss, 
particle re-acceleration processes and bubble dynamics in 
this radio source. The spectral slope of radiation emitted 
from different parts of the source is indicative of radiative 
ageing (or lack of it) of emitting particles. The low and high 
frequency spectral index maps are presented in Fig. [6] These 
maps show that between 240 and 621 MHz the spectral in- 
dex distribution is fairly uniform with no strong steepening 
or gradients across the source (a mean = —1.02 ± 0.17 for 
south-western bubble and a mean — —0.95 ± 0.10 for the 
north-eastern bubble, which are marked in Fig. [7}, implying 
a lack of strong radiative losses between these frequencies 



(Fig. |6j left panel). However, the high frequency spectral 
index map between 1.4 and 4.8 GHz presents a strikingly 
different picture (Fig.[5]right panel). On this map, while the 
south-eastern bubble still has the same average spectral in- 
dex value around -1 (« mean = —1.06 ± 0.15) - implying a 
straight synchrotron spectrum between low and high radio 
frequencies - the north-west bubble, in contrast has devel- 
oped an extremely steep spectrum (a mean = — 1.6 ± 0.20), 
most likely due to strong radiative energy losses in this part 
of the source. 

Such a situation might be explained if we assume that 
the north-western bubble showing the very steep spectrum 
is buoyantly rising away and then detaching itself from 
a centrally located mechanism of energy injection, while 
this source or mechanism is still active and possibly in- 
jecting relativistic particles into the relatively flatter spec- 
trum south-eastern bubble. One possibility is AGN related 
activity connected to a super-massive nuclear black hole, 
which is suggested by the large bolometric radio luminos- 
ity of MRC 0116+111 (L ra dio = 3.64 x 10 34 W, see Ta- 
ble 2), though we fail to detect a nuclear core or active 
jets (Fig. [1]). It is plausible that these bubbles are trac- 
ers of a previous cycle of AGN activity. Thus, a possible 
scenario is that these bubbles were inflated by the pres- 
sure of radio-jets in a previous very energetic episode of 
AGN activity which has either stopped now or has be- 
come too feeble to be detectable. During the later stages 
of evolution the dynamics of the remnants of the radio- 
jets will be domina t ed by buoyancy forces and drag (e.g. 
iGull and Northoverl (|l973T ). IChurazov et al l i|200lD ). as they 
interact with the cluster atmosphere in the inner core region. 
Our detection of a possibly buoyantly rising plasma bubble 
is also significant from theoretical point of view, as one of 
the most appealing and widely applied models in literature 
for heating of ICM involves buoyant gas bubbles, inflated 
by AGN source, t hat subsequently rise through the ICM 
and h e at it up (e.g. 



iChurazov et al 
il2002l); iBruggen and Kaiseij (1200: 



20021) ;lBohringer et all 
Bruggen et alj (I2002T) ; 



ICattaneoand Tevssied (120071 ); iRuszkowski et al.l l|2008h ) 
We further discuss the physics of these radio bubbles in the 
sections below. 



3.5 Remarkable morphological similarity with the 
twin radio-bubbles of M87 (Virgo A) 

Our radio observations presented above clearly show that 
MRC 0116+111 is a striking example of a radio-halo source 
consisting of twin plasma bubbles, possibly inflated in a 
previous episode of AGN outburst from a super-massive 
black hole. We wish to point out that, in this respect, 
it closely resembles the radio emitting bubbles observed 
around M87, the dominant central radio galaxy of the Virgo 
cluster (lOwen. Eilek and Kassirn] [ioOO) . 3C84, the central 
cD in Perseus cluster (jPedlarl Il990l) and 3C218 the cen- 
tral radio source in Hydra cluster ( McNamara et aI1l2000l ; 
I Wise et al.ll2007i ). However, as we have pointed out above, a 
notable major difference is that while these (and other clus- 
ter centre radio mini-haloes, see Table 3) exhibit the central 
galaxy's active nucleus (central core) and visible jet struc- 
ture, MRC 0116+111 in contrast is probably a quiescent 
source, lacking an active AGN and ongoing radio-jet activ- 
ity. Therefore, MRC 0116+111 is quite unusual and presents 
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Table 3. Comparison of MRC 0116+111 with other known mini 
radio haloes^. Column 1: Cluster name. Column 2: Cluster red- 
shift. Column 3: Log X-ray luminosity of the cluster in energy 
range 0.1-2.4 keV. Column 4: Log 1.4 GHz radio luminosity. Col- 
umn 5: Log radio halo radius in kpc. 



Name 


Redshift Log L x Log Li.4GHz 


Log -Rh 




(erg s" 1 ) (W Hz" 1 ) 


(kpc) 



Perseus 


0.018 


44.82 


24.27 


2.12 


A2390 


0.228 


45.13 


24.99 


2.26 


A2626 


0.060 


44.03 


23.36 


1.85 


RX J1347.5-1145 


0.451 


45.65 


25.28 


2.41 


Z7160 


0.258 


44.93 


24.34 


2.24 


RBS 797 


0.350 


45.31 


24.63 


2.01 


M87 


0.004 


43.90 


24.72 


1.60 


MRC 0116+111 


0.131 


< 43.15* 


24.66 


2.08 



Data from lCassano. Gitti and B runctti ( 2 0081) and lBirzan et al] 
20081 1 

* Upper limit estimated by us from ROSAT All Sky Survey 
(RASS) data in 0.5 - 2 keV band 




Figure 6. Low and high frequency spectral index maps of MRC 0116+111 obtained by combining the GMRT radio maps at 240 and 
621 MHz (left panel) and VLA maps at 1425 and 4860 MHz (right panel). Both pairs of maps have the matched resolution of 12" 
(FWHM). To prevent spurious structures from appearing, only pixels above ~3.5 times the rms noise level were included by giving cut-offs 
at values 4 mjy/beam, 0.5 mjy/beam, 0.27 mjy/beam and 0.15 mjy/beam at 240, 621, 1425 and 4860 MHz respectively. The colour-bar 
shows the spectral index values. A significant spectral steepening (mean a = —1.6) for the north-western bubble at high-frequency can 
be clearly seen. 



a unique opportunity to understand the dynamical evolution 
and impact of radio jet inflated plasma bubbles on the ICM, 
when the central AGN activity approaches cessation (prob- 
ably when the fuel supply to the SMBH is exhausted), or 
when the non thermal bubbles start rising buoyantly and 
detaching themselves from the flow. 

For a closer comparison, in Fig. [7] we show the 621 
MHz radio image of MRC 0116+111 and a deep VLA 



90 cm image of M87 published by lOwen. Eilek and Kassinil 
(2000) (the original image has been reflected and rotated 
for better comparison with MRC 0116+111). The morpho- 
logical similarity between the t wo radio sources is striking. 
lOwen, Eilek and Kasslml (|2000h found that in M87, beyond 
the well known inner-jet and lobe region (which appears as 
the 2 kpc scale red-orange patch near the centre of M87, 
Fig. , there is diffuse outer structure extending upto ~40 



10 Bagchi et al. 






Figure 7. Upper figure: GMRT 621 MHz false-colour im- 
age of MRC 0116+lll(beam: 6" FWHM circular, plotted in- 
side the box) showing the twin bubble-like diffuse radio struc- 
ture. The color-bar shows surface brightness in mjy/beam units 
and position of brightest central galaxy (cD) is shown by a 
symbol. The dotted lines delineate the two radio bubbles hav- 
ing projected dimensions ~27.5" X 17.5" (semi-major and minor 
axes of ellipse) for the smaller elliptical north-west bubble and 
~35" (circular radius) for the spherical larger bubble south-east 
of the cluster centre. Lower figure: The VLA 90 cm radio 
image of M87, the dominant central radio galaxy in the Virgo 
Cluster ijOwen. Eilek and Kassiml [20001. The original map has 
been reflected and rotated for a better comparison with the 
MRC 0116+111 radio map (reproduced by permission of the 
AAS). 



kpc from the nucleus, which is known to host a central super- 
massive black hole with a mass of 3.2 x 10 9 M Q ijHarms et all 
1994). Two large scale collimated flows emerge from the 
inner-jet region, one directed north-eastward and the other 
oppositely directed to south-westward (orientations refer to 
the rotated image). Particularly intriguing is a well-defined 
pair of edge-brightened, 'ear-shaped' torus-like structure at 
the end of south-west flow, which is reminiscent of a sub- 
sonic vortex ring. The north-eastern flow develops a gradual 
but well-defined S-shaped southward twist, starting only a 




Figure 8. An image processed view of GMRT 621 MHz image 
of MRC 0U6+U1 shown in Fig. [TJ We have used the 'Sobel' 
gradient filter operation as implemented in AIPS task NINER to 
enhance sharp intensity gradients and suppress diffuse emission. 
Highly intricate filamentary structures of radio emitting plasma 
interior to the radio-halo can be seen. 



few kiloparsecs beyond the inner lobes. Finally, both flows 
are embedded in a giant diffuse radio structure that might 
be described as two overlapping 'bubbles', each extending 
about 40 kpc from the nucleus. After reaching the halo, 
the flows gradually disperse, particularly the north-eastward 
flow, and appear to be filling the entire halo with radio-loud, 
filamented plasma. 

We observe similar flow pattern in the central region 
of MRC 0116+111, despite the limited spatial resolution 
available. As shown in Fig. [7] an edge brightened torus-like 
'mushroom' structure is observed about 40 kpc west of cen- 
tre of the larger south-eastern bubble, which might well be 
analogous to the peculiar 'ear-shaped' vortex structure of 
M87. Here the flow pattern sharply turns northward and 
appear to be flowing into the smaller radio bubble to the 
north-west. We observe an S-shaped flow pattern to north- 
east of centre, which further bends to south, and clearly re- 
sembles the filamentary structure found in the southern bub- 
ble of M87. In order to better visualize a possible filamentary 
structure of synchrotron plasma bubbles, we have digitally 
image processed the 621 MHz map of MRC 0116+111 with 
the AIPS task NINER for a gradient filter operation, which 
enhances sharp intensity gradients and suppresses diffuse 
emission. The resulting image reveals (Fig.[S|, akin to M87, 
a highly complex filamentary structure of plasma embedded 
in the radio-halo. However , with the available angular res- 
olution (6" FWMH beam size) far less details are visible in 
MRC 0116+111 as compared to M87, mainly because it is 
much farther away from us than M87. Of course, as men- 
tioned above, in our radio observations of MRC 0116+111 
we find no evidence for a strong AGN core, inner-jets and 
lobe structures, which clearly are present in M87. One fea- 
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ture which makes the case of MRC 0116+lllespecially in- 
teresting is that the inner plasma outflow structures in both 
these radio sources are surrounded by a pair of diffuse radio 
emitting 'bubble' shaped giant lobes. We note that large mis- 
alignment seen in these two cases, between the inner outflow 
structure and the surrounding mini-halo may in fact be fairly 
common, exhibited by several other mini-haloes. Examples 
include A2626 and the protot ypical mini-halo in Perseus-A 
l)Gitti et al.l2004l : |Pedlarill99fi l . All this hints towards a com- 
mon mechanism, namely the radio-jets from AGN, strongly 
interacting with the confining environment of a cooling core, 
forming an amorphous mini-halo structure. 

Moreover, it is also interesting to note that although 
MRC 0116+111 closely resembles M87 in its radio mor- 
phology, at ~240 kpc its total linear size is about three 
times larger than M87, which sp ans a scale of just ~80 kpc 
|Owen. Eilek and Kassiml (|200Ch and Table 2). The bolo- 
metric radio luminosity of both is about the same; I/radio = 
3.64 x 10 34 W for MRC 0116+ 111 and L radio = 3.6 x 10 34 
W for M87 l|Birzan et alj|2008l l. This in turn implies that 
the average volume emissivity of MRC 0116+111 is almost 
30 times smaller than that of M87. 



4 DISCUSSION 

The similarity of large scale radio plasma outflow structures 
observed in MRC 0116+111 and M87 suggest that both 
sources might have originated in, and their evolution gov- 
erned by, similar physical process and conditions prevailing 
in the central regions of their host clusters. The origin of 
extremely complex diffuse radio structures observed in the 
outer radio bubbles of M87 and other similar radio mini- 
halo objects is n ot yet fully understoo d. From Hydrody- 
namic simulation IChurazov et all (2001) suggested that the 
twin bubbles in M87 are buoyant bubbles of cosmic rays and 
magnetic fields, inflated by jets launched during an earlier 
nuclear active phase of the central galaxy, which rise through 
the co oling gas at roughly half i t s soun d speed. On the con- 
trary, lOwen. Eilek and Kassiml (|2000T ) favor a scenario in 
which the radio halo of M87 is not simply a relic of previous 
episode of AGN activity but presently it is 'alive', such that 
it is being supplied with relativistic particles coming from 
the central AGN and the innermost radio jets. As revealed in 
our radio observations, non-detection of a strong AGN core, 
inner-jets and lobe/hotspot structure in MRC 0116+111, 
coupled with the detection of a ~100 kpc scale radio bubble 
of very steep spectrum north-west of the cl uster centre sup- 
ports the 'buoyantly-rising bubble' model of lChurazov et all 
|200lT ). The flattened 'mushroom' shape of this plasma bub- 
ble indeed resembles a rapidly rising vortex-ring (Fig. 0), 
into which an initially spherical bubbl e will naturally trans- 
form due to viscosit y and drag forces (IChurazov et al.ll200ll ; 
Brug gen et al.| [2002) . The present quiescent /low state of the 
central engine suggests an on-off activity cycle of the su- 
per massive black hole. That such activity is intermittent 
and the central AGN is periodically rejuvenated is sug- 
gested by the high incidence of bubbles found in cluster 
cores that require heatin g to counter balance the cooling flow 
ijDunn and Fabi an 2006), and by the train of ghost-bubbles 
and pressure-ripp l es observed in Perseus and Hydra clusters 
jWise et al.ll2007l ; ISanders and Fabianll2007h . Also, the dis- 



tant cooling core cluster RBS797 shows an inner radio-jet 
pair oriented perpendicular to that of the extended radio- 
halo structure and X-ray cavities, suggesting re-s tarting ac- 
tivity of AGN (|Gitti. Feretti and Schindleij|2006T) . 

An important question raised by our observations is: 
what energizes the relativistic particles in the central re- 
gion of MRC 0116+111 which apparently counter balances 
the strong radiation losses and prevents severe steepening of 
radio spectrum? The properties and physical processes gov- 
erning the origin and evolution of radio mini-halo sources 
are poorly understood due to their rather small numbers 
known till now (les s than 10). Based on a sample of 5-6 mini 
haloes, iGitti et ail i|2004l ) have reported a positive correla- 
tion between the radio power of mini-halo and the cooling 
flow radiative power. From this, and because of their esti- 
mate of the lifetimes of the radiating electrons falling con- 
siderably short of the diffusion time-scale, they have argued 
that in mini-haloes the cooling flow (through the compres- 
sional work done on the ICM and amplification of frozen-in 
magnetic field) energizes the particle acceleration process 
through magneto-plasma turbulence acting on relic electron 
population probably left behind by the past episodes of AGN 
activity. This suggests a direct connection between cooling 
flows and radio mini-haloes. With a detailed spectral fitting 
of the radio spectrum and emission profile with the param- 
eters derived from the turbulent re-acceleration model, it 
will be possible to test the viability of this mechanism for 
the present source. 

An alternative proposal is that the radiating elec- 
trons in the mini-haloes are of secondary hadronic ori- 
gin (via the pion de cay route), created in proto n-proton 
(PcR-Pth) collisions (|Pfrommer and Ensslinl l2~004h of high 
energy cosmic ray (CR) protons (Pcr) with target ther- 
mal protons (Pth)- The CR protons may originate in sev- 
eral ways: in large-scale struct u re formation shock waves 
JColafrancesco and Blasi 1 19981 : lLoeb and Waxmanl [2000; 
iRvu et al.1 120031 : iBagchi et alj 120061 ), or thro ugh super- 
novae and galactic winds from central galaxies (Volk et al. 
Il99fj). or via powerful AGNs injecting CRs into the ICM 
Enss lin et al.lll997l : lAharonian 2002; P frommer and Ensslinl 



20 



2004). A clear signature of this process is emission of 7-ray 
photons of broad energy range ~ 100 MeV - 10 TeV, ac- 
companying decay of neutral pions (tv°) produced in inelas- 
tic proton-proton collisions. Therefore, the hadronic model 
should be testable in near future with sensitive 7-ray o bser- 
vations. Recently, iHinton. Domainko and Pope! (|2007n ap- 
plied the hadronic CR injection model to three galaxy clus- 
ters (Hydra-A, MS 0735.6+7421 and Hercules-A) harboring 
the most energetic AGN outbursts known (with mechani- 
cal energy input of AGN at least 10 61 erg). They concluded 
that the 7-ray observability of these objects is at the thresh- 
old of detectability with current (FERMI, HESS), and near 
future 7-ray instruments, if the observed radio bubbles are 
dominated by the pressure of relativistic particles. There- 
fore, such observations are obviously of great importance 
for deciphering the connection between the energetic CR 
and thermal energy content of clusters harboring the radio- 
haloes. 

We have estimated the volume of the two radio bubbles 
from their projected dimensions on the GMRT 621 MHz 
image (Fig. [7] Table 2). The volumes are: Vi = 4/3(7rr?) = 
6.78 x 10™ cm 3 for the spherical south-west bubble of ra- 
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dius n = 82 kpc, and V 2 = 4/3(7ra 2 fe) = 2.08 x 10™ cm 3 for 
the ellipsoidal north-west bubble (we model it as an oblate 
spheroid with semi-major axis a = 64.4 kpc and semi-minor 
axis b = 41 kpc). The volume V = Vi + V 2 = 8.86 x 10 70 
cm 3 is the total volume occupied by relativistic particles and 
magnetic fields, which we assume to be mainly supplied by 
the radio jets and not produced in local shocks. The linear 
dimension of the present radio bubbles is hubble ~ 100 kpc, 
which is much larger compared to the typical dime n sions of 
X-ray cavities (Icavity ~ 1 - 20 kpc; iBirzan et all Hooi)), 
and they are comparable to the ~100 - 200 kpc size su- 
per cavities found in MS0735 . 6+7421 and Hydr a clusters 
i|McNamara et all (|2000l , 120051 '); I Wise et all (|2007h V 

The large radio luminosity and linear size of the radio 
bubbles in MRC 0116+111 suggest that radio jets from AGN 
must have done large amount of mechanical work (pdV) in 
inflating them against the external pressure of thermal ICM 
gas, as well as injected enormous non-thermal plasma en- 
ergy (magnetic field and relativistic particles), which has 
remained stored inside the bubbles, as they rose buoyantly 
under pressure equilibrium with the external medium. Fu- 
ture X-ray observations of MRC 0116+111 would reveal if 
the relativistic fluid filling these bubbles has pushed aside 
the thermal gas, excavating giant cavities in the X-ray emit- 
ting ICM. We have detected a buoyantly rising radio bub- 
ble ~ 100 kpc from the centre of radio-halo (Fig. |S} . Such 
bubbles result from displacement of the thermal gas, which 
creates a low-density bubble maintained in pressure bal- 
ance with the surrounding medium and rising subsonically. 
These non-thermal plasma bubbles associated with X-ray 
cavities found in the ICM of clusters are turning out to 
be excellent 'calorimeters' for guaging the total mechani- 
cal power of radio-jets from AGNs. Such measurements are 
extremely useful as they not only preserve the past record 
of enormous outbursts of the central SMBH, but are also 
independent of the radio properties of the jets (i.e., their 
present radiative output and magnetic field strength) , allow- 
ing a good estimate of energy inject ed into the surroundin g 
thermal plasma for heati ng it (e.g.. [Churazov et al.l ((2002); 
IBirzan et all (|2004 120081 ')'). 

A thermodynamical measure of the total energy re- 
quired to slowly inflate a hot cavity /bubble in ICM is given 
by its enthalpy, H ctrt i=[7/(7 — l)]pV, where 7 is ratio of 
specific heats of the gas inside the cavity, p is the pres- 
sure of_e2rte£rial_jnedlum £ind_y_istlie volume of the cav- 
ity l|Churazov et all |2002| ; IBirzan et all 120041 '). For a radio 
jet inflated non-thermal bubble, which has completely dis- 
placed the thermal medium, 7 = 4/3 and H cav = 4pV. The 
external pressure is p — picukT '/(jjms) (where picM is 
the density of ICM, T is the gas temperature, fi—0.6 the 
mean molecular weight, and tuh is mass of the hydrogen 
atom), which is presently not known for MRC 0116+111. 
We can obtain an order of magnitude estimate for H cav if 
we take fiducial values, picm = 1-67 x 10 -26 gm cm -3 (for 
a proton number density of 10 -2 cm -3 ) and T — 2.32 x 10 7 
K (2 keV), which are typical of the gaseous environment 
around a mini radio-halo sour ce like M87 or Perseu s-A at 
the centre of a cooling core (|Ghizzardi et al.l I2004F ) . This 
gives a pressure p = 5.3 x 10~ n dynes cm -2 and an en- 
thalpy H C av = 2.06 x 10 60 V7o erg for a cavity volume 
V70 = 17(70 cm 3 ). 

Using the measured volumes of the two radio bubbles 



in MRC 0116+111 (see above), we estimate an enthalpy 
Hi = 1.39 x 10 61 erg for the south - east bubble, and 
H2 — 4.28 x 10 60 erg for the north - west 'buoyant' bub- 
ble, resulting in total enthalpy -H1+2 = 1.82 x 10 61 erg. 
We require X-ray data to estimate the external ICM pres- 
sure around the bubbles; but conservatively, even if this 
is 1/10 of the value assumed above, we get a total en- 
thalpy i/i+2 ~ 10 60 erg which is still quite large, but for 
an order-of-magnitude consistent with the large radio lu- 
minosity and huge volumes of these bubbles. The energy 
injection rate (mechanical luminosity) of radio jets, that 
inflated the radio bubbles, is Lj « Hi+2/tj, where tj is 
the time interval for which jets were active out of the to- 
tal activity cycle time td of the AGN. Various lines of evi- 
dence show that for radio galaxies typically, td «fewxl0 8 
y, and assuming tj ~fewxl0 7 y wit h J/1+2 ~ 10 60 erg 
would lead to Lj -TewxlO 44 erg s' 1 . IBirzan et all (|2008T ) 
derived an empirical relation (with considerable scatter); 
log = (0.35 ± 0.07) log Liagh* + (1-85 ± 0.17), con- 
necting the estimated mechanical jet luminosities of radio 
AGNs {Lj = L cav ) that have associated X-ray cavities, and 
their 1.4 GHz monochromatic radio luminosities. In this re- 
lation Liaghz is in units of 10 24 W Hz -1 and the cavity 
(jet) luminosity L cav is expressed in units of 10 42 erg s _1 . 
From this relation, and with Li.4GHz = 4.57 x 10 24 W Hz" 1 
for MRC 0116+111 (Table 2), we obtain Lj = 1.2 x 10 44 
erg s _1 , which is consistent with the earlier 'calorimetric' 
estimate. Moreover, adopting tj = 10 7 — 10 8 y, the same as 
the synchrotron age derived in Sect. 3.3, we estimate that 
the total energy content (enthalpy) of the bubble sytem is 
~4x (IO 08 — 10 59 ) erg, which is less than our earlier esti- 
mate (10 60 - 10 61 ) erg, but we believe it represents best 
the situation in MRC 0116+111. However, the exact figure 
can only be obtained from the future X-ray observations of 
MRC 0116+111. 

Indeed, both observations and theory suggest that the 
energy input from AGN into the cluster atmosphere at 
this rate is sufficient to heat i t up and balance the coo l- 
ing loss l|Raffertv et all 120061 ; iDalla Vecchia et all 120041 '). 
This may even halt (or reduce drastically) the accretion 
flow of matter onto t he SM BH in some systems (e.g., 
ICattaneo and Tevssierl |2007)) - thus cutting off the fuel 
supply to the 'central engine' and completing the negative 
feedback loop. We wonder, if this AGN mediated feedback 
process could also be at work in MRC 0116+111, giving rise 
to its relaxed, quiescent looking appearance? If confirmed 
with future X-ray observations, such large enthalpy of the 
bubble system in MRC 0116+111 would make it harbor one 
of the most energetic radio outbursts known, comparable to 
the energetics o f the super cavities in MS0735. 6+7421 a nd 
Hydra clusters (|McNamara et alJliooS IWise et alj|2007f ). 

The ultra steep spectrum north - west bubble is located 
at a projected distance of R~100 kpc from the central cD 
galaxy, which is inferred to be buoyantly rising in the ICM. 
We do not know how far back in time it was inflated and 
where was it injected into the ICM. Assuming that it was ini- 
tially inflated rapidly by an AGN-fed radio jet near the clus- 
ter centre, and subsequently it detached and rose subsoni- 
cally at ~ 1/2 the speed of sound C s = 717.6 [T/(2 keV)] 1/2 
km s~\ it will take t s/2 = 2C' 1 R = 2.7 x 10 8 y for it to 
arrive at its present location (neglecting projection effects). 
This time-scale is much longer than the electron cooling time 
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scale t cool w 2 x 10 7 y (taking B = 10 fiG] at 4.86 GHz, 
thus consistent with its observed very steep radio spectrum 
(Fig. [6]). However the actual rise time should be somewhat 
shorter than the above estimate when we take into account 
the initial rapid inflation and acceleration phase (due to for- 
ward momentum of the radio jet), at the end of which the 
bubble will attain a constant terminal speed Vterm and rise 
in pressure balance with the external medium. Another rel- 
evant time-scale is the 'buoyant rise' time tbuoy, which is 
the travel time for a bubble moving with the terminal speed 

Vie 



The X-ray data show that u sually tbuoy ~ 2 x i. 



/2 



ijBirzan et al.ll2004l : Iwise et al.ll2007h 

Several galaxy clusters also show X-ray surface bright- 
ness depressions that have no obvious association with 
bright radio emission - the s o called 'ghost cavities' , such 
as those seen in Abell 2597 (jMcNamara et al.ll200ll ). and 
the matched pair of cavities in Hyd ra and Perseus clusters 
( Fabi an et al l200d : IWise et al.ll2007h . These depressions are 
thought to have been created by AGN outbursts that oc- 
curred in the more distant past, but whose radio emission 
has faded over time due to rad i ation and expansion losses 
(eg. lEnsslin and Gopal-Krishnal (|200ll l). An interesting pos- 
sibility, suggested by our discovery of an ultra steep spec- 
trum detaching radio bubble in MRC 0116+111 is that, that 
such buoyantly rising radio bubbles are precursors of the 
ghost cavities observed in many X-ray observations. Such 
ghost bubbles provide important clues for understanding 
how AGNs heat cluster atmospheres. A crucial issue is the 
survival of detached radio bubbles upto large distances from 
the central source. The stability and dynamics of radio bub- 
bles, once the jets have turned off (as in the present radio- 
halo MRC 0116+111) is a topic of great importance for the 
interpretation of X-ray observations of clusters. Chandra ob- 
servations show that some of the ghost cavities are able to 
reach ~20 - 100 kpc distances from the central AGN, while 
it is found in several numerical simulations that hot ther- 
mal bubbles, once released into the cluster atmosphere at 
the end of a terminating radio-jet, initially rise supersoni- 
cally but are rapidly fragmented by Ray leigh- Taylor (RT) 
and Kelvin-Helmholtz (KH) instabilities. To prevent their 
destruction, the role of magneti c field in st a biliza tion of the 
rising bubbles was e xplored by |Pe Yound (|2003T) . Numeri- 
cal simulations (e.g., iRuszkowski et aTP 20081 )) find that in 
an upward-rising cosmic-ray filled bubble, the magnetic field 
'draping' effect at the upper contact surface of bubble indeed 
has a strong stabilizing effect which prevents cross-field dif- 
fusion of particles. A substantial fraction of cosmic-rays can 
thus be confined inside the bubbles on buoyant rise time- 
scales (~ 10 7 " 8 y) even when the parallel diffusivity coeffi- 
cient is very large. Magnetic field stabilization of the radio 
bubbles in MRC 0116+111 is an interesting possibility worth 
exploring, as these are still intact quite far away from the 
central cD: ~100 kpc for north - west bubble and ~40 kpc 
for south - east bubble, requiring some stabilizing mech- 
anism for their survival. Detailed polarization mapping of 
these features is required to check if the theoretical 'draped' 
structure of protecting magnetic field does exist. 

Lastly, MRC 0116+111, with its dual bubble-like mini 
radio-halo morphology, provides an excellent test-bed for 
understanding the dynamics of radio-jet inflated buoyant 
plasma bubbles in the hot cluster atmospheres, and the 
role of AGN feedback in transporting and mixing of met- 



als in the ICM. Low-entropy, metal enriched cold gas may 
be uplifted from the cluster centre along with the energetic 
wakes of AGN-driven plasma bubbles upto ~ 10 — 100 kpc 
scales, which undergoes heating and mixing with ICM via 
hydro dynamical turbulen ce. Recently, from HST observa- 
tions iFabian et al.l (|2008T ) discovered an extensive network 
of H a emitting cold, apparently magnetically supported, fil- 
aments surrounding the mini radio-halo Perseus-A. They in- 
ferred that these cold filaments are dredged up from the cen- 
tre of the galaxy by the radio-emitting bubbles, buoyantly 
rising in the hot intra-cluster gas, and act as markers of the 
AGN feedback process. 

We expect that future more detailed radio, X-ray, 7- 
ray and optical emission line studies of the galaxy cluster 
hosting MRC 0116+111 would provide crucial data, that 
would allow us to gain a much better understanding of the 
complex thermal and non-thermal plasma processes taking 
place in the gaseous environment around this puzzling yet 
extremely interesting radio-halo source. 



5 SUMMARY AND CONCLUSIONS 

(i) We have presented detailed radio and optical observa- 
tions of MRC 0116+111, which clearly have revealed a lumi- 
nous, steep-spectrum, diffuse radio-halo source of ~ 240 kpc 
dimension located at the centre of a X-ray weak, low-mass 
cluster of galaxies at redshift z = 0.131. 

(ii) On larger scales (~100 kpc) the radio-halo emission 
comprises of two diffuse 'lobe'-like structures surrounding 
the central cD and its associated group of galaxies, which 
we have interpreted as being a pair of radio emitting gi- 
ant plasma bubbles, filled with relativistic particles (cosmic 
rays) and magnetic field. 

(iii) We do not detect any ongoing AGN activity, such 
as a compact core or active radio jets feeding the plasma 
bubbles. Thus, it is possible that these bubbles are tracers 
of a previous cycle of AGN activity and were inflated by the 
pressure of radio-jets in a highly energetic episode of AGN 
outburst, which has either ceased by now, or has become 
too feeble to be detectable. 

(iv) The physical size, radio spectral index and morphol- 
ogy of MRC 0116+111 are broadly consistent with the prop- 
erties of a cluster radio mini-halo, which are rare objects, 
have low surface brightness and possess a steep spectral in- 
dex (a<^-l). They are mostly observed at the centres of 
cooling-core clusters, which indicates that their origin and 
evolution are closely linked to the energy feedback from the 
central AGN via radio jets into ICM and its cooling/heating 
processes. The radio luminosity of MRC 0116+111 at 1.4 
GHz, 4.57 x 10 24 W Hz" 1 , and the bolometric radio lumi- 
nosity (over 10 MHz - 10 GHz range), 3.64 x 10 34 W, would 
place it amongst the most luminous radio haloes known. 

(v) The integrated radio spectrum shows a sudden steep- 
ening beyond the 'break' frequency ~400 MHz. The low- 
frequency spectral index is — 0.55(±0.05) and the high- 
frequency spectral index is — 1.35(±0.06). The electron spec- 
tral age was derived to be (1.33 — 0.64) x 10 s y, for magnetic 
field values 1-10 piG. 

(vi) The pair of huge (~100 kpc scale) plasma bubbles 
show marked spectral index differences. Between 240 and 
621 MHz (GMRT) both bubbles show a similar mean spec- 
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tral index a mean « —1, indicating a lack of strong ra- 
diative losses. However, whereas between 1.4 and 4.8 GHz 
(VLA) the south-eastern bubble still has the same mean 
spectral index of — 1.06 ± 0.15, the north-west 

bubble, in contrast, has developed an extremely steep spec- 
trum ctmean = —1.6 ± 0.20, clearly due to strong radiative 
energy losses. We argued that ongoing in-situ particle re- 
acceleration, probably via merger induced shocks and/or 
turbulence is taking place in the south-eastern bubble, and 
the north-western bubble is detaching and buoyantly rising 
away from the centre in the hot cluster atmosphere. Such up- 
rising, ageing bubbles are potential precursors of the giant 
X-ray dark cavities and could inject ~ io 60-62 erg mechani- 
cal energy into the ICM for heating it, as observed in X-ray 
observations of several clusters. 

(vii) We discovered a remarkable morphological similar- 
ity between MRC 0116+111 and M87, the well known dom- 
inant central radio galaxy in the Virgo Cluster. Similar in- 
ternal plasma outflow structures and network of magnetic 
filaments surrounded by a pair of giant outer radio bubbles 
are observed. This suggests that both these sources might 
have originated in, and their evolution governed by, similar 
physical processes and conditions prevailing in the central 
regions of their host clusters. Inspite of the notable similari- 
ties, at ~240 kpc, the linear size of MRC 0116+111 is about 
three times larger than M87. The bolometric radio luminos- 
ity of both is about same: 3.6 x 10 34 W. This in turn implies 
that the average volume emissivity of MRC 0116+111 is 
roughly a factor of 30 smaller than that of M87 halo. 

(viii) From the known empirical correlation between the 
estimated mechanical jet luminosities of radio AGNs that 
have associated X-ray cavities, and their 1.4 GHz monochro- 
matic radio luminosities, we estimated the total enthalpy 
(free energy) of bubbles in MRC 0116+111, ~ 10 59 - 60 erg 
and the mechanical luminosity of radio jets which inflated 
them, to be ~ 10 44 erg s _1 . Energy input from a central 
AGN into the cluster atmosphere at this rate is sufficient to 
heat it up, drive a massive outflow from the AGN, which may 
counter balance a possible cooling loss, and may even halt 
(or reduce drastically) the accretion flow of matter onto the 
supemassive black hole, thus cutting off the fuel supply pow- 
ering the central engine and completing the feedback loop. 
This AGN mediated feedback process in MRC 0116+111, 
could explain its relaxed, quiescent looking appearance. If 
confirmed with future X-ray observations, such large en- 
thalpy of the bubble system in MRC 0116+111 and its possi- 
ble AGN-feedback cycle would make it host one of the most 
energetic radio outbursts known, comparable to the ener- 
getics of the super cavities in MS0735. 6+7421 and Hydra 
clusters. 

(ix) MRC 0116+111, with its rare, twin bubble-like mor- 
phology of mini radio halo and possessing other striking 
properties as revealed in our observations, provides an ex- 
cellent opportunity to understand the energetics and the 
dynamical evolution of radio-jet inflated plasma bubbles in 
the hot cluster atmosphere. It also help in understanding 
what possible role AGN-feedback plays in heating of ICM, 
transporting and mixing of metals in the ICM, the evolution 
of central galaxies and the growth of nuclear supermassive 
black holes. 
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